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Abstract. 

■^1- ' Exact solution is obtained for electromagnetic field around a conducting cylinder 

of infinite length and finite radius, with a periodical axial current, when the wave 
length is much larger than the radius of the cylinder. The solution describes simul- 

O^l ' taneously the fields in the near zone close to the cylinder, and transition to the wave 

(-H , zone. Proper long- wave oscillations of such cylinder are studied. The electromagnetic 

energy flux from the cylinder is calculated. These solutions could be applied for 

description of the electromagnetic field around relativistic jets from active galactic 

fvj nuclei and quasars and particle acceleration inside jets. 
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1. Introduction 

q ■ Objects of different scale and nature in the universe: from young and 

very old stars to active galactic nuclei (AGN) (Bridle, 1984),(Eilek, 
1984),(Eilek and Hughes, 1990) show existence of collimated outbursts 



S^. • - jets. Geometrical sizes of jets lay between parsecs and megaparsecs. 

The origin of jets is not well understood and only several qualitative 
mechanisms are proposed. Theory of jets should answer to the question 

-*— > ■ of the origin of relativistic particles in outbursts from AGN, where 

synchrotron emission is observed. Relativistic particles, ejected from 
the central machine rapidly loose their energy so the problem arises 
of particle acceleration inside the jet, see reviews (Begelman et al., 

X ' 1984), (Bisnovatyi-Kogan, 1993). 

C^ ■ It is convenient sometimes to investigate jets in a simple model 

of infinitely long circular cylinder (Chandrasekhar and Fermi, 1953). 
Magnetic fields in collimated jets is determining its direction, and axial 
current is stabilizing its elongated form at large distances from the 
source (AGN) (Bisnovatyi-Kogan et al., 1969). When observed with 
high angular resolution these jets show structure with bright knots sep- 
arated by relatively dark regions (Bridle, 1984), (Eilek, 1984), citetmw. 
High percentages of polarization, sometimes more then 50% in some 
objects, indicates the nonthermal nature of the radiation which is well 
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explained as the synchrotron radiation of the relativistic electrons in a 
weak but ordered magnetic field. Estimation of the life time of these 
electrons, based on the observed luminosities and spectra, often gives 
values much less the their kinematic ages tk/c, where d is the distance 
of the emitting point from the central source. Because the jet flow most 
likely originated from an outburst or continuous outflow from the cen- 
tral source, there is a necessity of continuous reacceleration of electrons 
in the jets in order to explain the observations. Acceleration mechanism 
for electrons in extragalactic jets, proposed in (Bisnovatyi-Kogan and 
Lovelace, 1995) considers that intense long-wavelength electromagnetic 
oscillations accompany a relativistic jet as a result of the non-steady 
mechanism of the jet's generation in the nucleus of the source. The 
electromagnetic wave amplitudes envisioned are sufficient to give in 
situ acceleration of electrons to the very high energies observed > 10 13 
cV. It was assumed that jets are formed by a sequence of outbursts 
from the nucleus with considerable charge separation at the moment of 
the outburst (Bisnovatyi-Kogan et al., 1969). The direction of motion 
the outburst is determined by the large-scale magnetic field. The out- 
burst are accompanied by an intense electromagnetic disturbance which 
propagates outward moving with the jet material in the direction of the 
large scale magnetic field. It was suggested in (Bisnovatyi-Kogan et al., 
1969) that a toroidal magnetic field, generated during the outbursts is 
important for the lateral confinement of of the jet. 

When the plasma density in surrounding medium is small, the elec- 
tromagnetic wave generated by the nonpotential plasma oscillations 
of the confined body is emitted outside and can accelerate particles. 
When the emitted wave is strong enough it washed out the medium 
around and the density can become very small, consisting only of the 
accelerated particles. The action of the oscillating knot is similar to the 
action of the pulsar, as inclined magnetic rotator. Both emit strong elec- 
tromagnetic waves, which could effectively accelerate particles (Pacini, 
1967), (Gunn and Ostriker, 1970). Long-periodic proper oscillations 
in the plasma cylinder with a finite radius, and emission of electro- 
magnetic waves had been studied in (Bisnovatyi-Kogan and Lovelace, 
1995), and in a simpler model in (Bisnovatyi-Kogan, 1996). Enhanced 
oscillations in such cylinder have been studied in (Bisnovatyi-Kogan, 
2004). Both models are represented below. 



2. Cylinder with oscillating current 

Consider infinitely conducting circular cylinder in vacuum. This model 
is valid at low density in surrounding plasma, which cannot screen the 
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emitting electromagnetic wave. Maxwell equations are (Landau and 
Lifshits, 1982) 

1<9E Air 
divB = 0, rotB = -^- + — j, (1) 

c ot c 

1 BT-i 

rotE = — , divE = 47rp e . (2) 

c ot 

For periodic oscillations with all values ~ exp(— iuot) they read 



divB = 0, rotB = -— E + — j, rotE = — B, divE = 0. (3) 

c c c 

We use the same definitions for all complex values depending on coordi- 
nates. Consider infinitely long cylinder with zero charge density, where 
in the cylinder coordinate system (r, (f>, z) the only nonzero components 
are E z , B^, j z , and d/d(p = d/dz = 0. Only two valid equations remain 
from the system (3): 

dE z ioj „ 1 d{rB&) iuj _ Air 

—^ + —B^ = 0, -^-r^ + —E z j z = 0; 4 

dr c r dr c c 

from which, we obtain equation for E z : 

\ d ( dE z \ J 2 ,_, 4niuj , . 

r dr \ dr I c z & 



3. Vacuum solution 

In vacuum j 2 = 0. Using non-dimensional variable x = ruj/c, we obtain 
from (4), (5) 

x 2 E' z + xE' z + x 2 E z = 0, 5^ = i££. (6) 

Here ' denote differentiation over x. The equation (6) belongs to Bessel 
type and has a solution 

E z = CiJo(x) + C 2 y (x), ^ = -i[d Ji(x) + C 2 yi(x)]. (7) 
Relations for Bessel functions are (Gradshtein and Ryzhik, 1964) 

JJ(x) = -Ji(x), 5tf(z) = -H(x). (8) 

General solution for physical values, with account of the time depen- 
dence, is obtained from the real part of the complex solution at 
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exp {-iujt) = cos ut-i smut, C x = c[ r) +iC[ i} , C 2 = C^ ] +iCf ] . 

(9) 
The general solution in vacuum is 

E z = [C[ r) j (x) + C| r) y (x)] cos cut + [Cf } J (x) + cJVo(x)] sinwt, 

(10) 

£ = -[C{ r) Ji(x) + cf Vi(x)] sinwi + [Cf } Ji(x) + cf } Yi(x)] coswi. 

(11) 
The boundary condition far from the cylinder follows from the demand 

that there exist only expanding wave. It means, that only functions 

depending on the combination (x — uot) survive. Using asymptotic of 

Bessel functions at large arguments (Gradshtein and Ryzhik, 1964) 

t i \ ~ [^ i 7^ ^ t ( \ ~ Z^" ■ i 7^ ^ 
Jo(x) ss \ — cos(x — -), J\\x) ~ W — 5111(1 — -), 

\ TTX 4 V 1TX 4 

/ 2 7T / 2 7T 

Yo(^) ~ V — sin(x — -), Fi(x) ~ — W — cos(x — -), at x>l. 

V 7TX 4 V 7TX 4 

(12) 

_ r {r) Ml) _ r (r) , 
Ly 2 i u 2 " u l ) ied 

to the following solution at large distances 



we obtain for the expanding wave C} = — C^ , C^ = C{ , leading 



£ 2 « J — [cf r) cos(x- j-wt) + cf ) sin(x- j-ut)], B (j> = -E z . 

(13) 
The general vacuum solution, satisfying conditions at infinity reads as 



E z = [cf r) J (x) + C^ r) Y (x)] cosujt + [-C$ ] J (x) + c[ r) Y (x)] smut, 

(14) 

5^ = -[C 1 (r) Ji(x) + cf ) yi(x)]sinwt+[-Cf ) Ji(x) + C 1 (r) yi(x)]coswt 

(15) 



4. Solution inside the cylinder 

The equation in the matter are 

x 2 E" + xE' + x 2 £ 2 + — x 2 j z = 0, 5a = iE' (16) 

LO 
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A solution of the non-uniform linear equation (16) is a sum of a general 
solution of the uniform equation, and a particular solution of the non- 
uniform one £q(x). 

E z = SiMx) + £ 2 Y (x) + £ Q (x). (17) 

The function Yq(x) is singular at x = 0, so for a finite solution £2 = 0. 
We look for a particular solution in the form E z = £{x)Jq{x). First 
order equation with respect to F = £' follows from (16) 

x 2 (F'J + 2FJ ) + xFJ + — x 2 j z = 0. (18) 

From this equation we obtain the general solution for the amplitude of 
the electric field in the matter, in presence of periodic EEF: 



E g = - — J (x)l -727-X / zJo(z)j g (z)dz + S 1 J (x). (19) 



4m [* dy v 

u Jo yJoiv) Jo 

Consider waves much longer than the radius of the cylinder tq 



x = — < 1. (20) 

c 

Than use expansion at x <C 1 (Gradshtein and Ryzhik, 1964), 

Jb«l-X' J i~l' Ibw-Inf, *«- — . (21) 

4 2 7T 2 7T2; 

Using (21) we obtain from (19) the solution for long waves 

_ 2iuj f x r , .dy , . 

£,= — 2- / 4 y — + ^1, 22 

where I 2 (r) = I z {y) is the complex amplitude of the electrical current 
inside a cylindrical radius r = cy/u 

rr „2 /^x 

I z = 2tt j z rdr = 2ir^ j z xdx (23) 

JO oj z Jo 

Complex values: the function I z {r) and the constant £\ are 

J, = lW + ilW fc^+tfj . (24) 
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5. Matching of solutions and long-wave limit 

The total electrical current through the cylinder Iq = I z {tq), and fields 
on its surface (inside) Eq = E z (tq), Bq = B^tq) (real parts of complex 
relations) are defined as 



(r) (i) 

Iq = Iq cos ut + Iq sin ut, Eq 



c z Jo y 



cos ut + 



c 1 h y 



dy_ 

V 



sinwi, 

(25) 



B, 



o 



r(r) 



2^Ij_ XQ p (i) 

C 2 x + 2 &1 



cos ut + 



2u J« 



X Mr) 



c 2 x 



£ 



sin uot. 



At small xq the external solution on the surface of the cylinder is written 

as 

Ez0 = [C{ r) + Ci r) - In ^1 cos ut + [-Ci r) + C[ r) - In ^1 sin ut, (26) 

7T 2 7T 2 



B f! 



i(r) Xq 



■Ar) 



(r)XQ 



i{r) 



cK r,~u _ c ^_] s[nujt + [_ C ^^i - C K {>—\ cos Ut 



2 ttxq 2 ~ vrxo' 

All field components are continuous at the cylinder surface in absence 
of the surface charges and currents. Matching magnetic and electrical 
fields we obtain using (25), (26) coefficients in the solution of the 
external electromagnetic field, which are determined by the periodic 
electrical current in the cylinder: 



\itxq vr 2 



2ul, 



(r) 



XqU 



C 2 Xq 



r"ir>(vA (27) 

Jo y 



a 



(r) 



r(') 



7TXQ 



^ ]u ^\ = ^Io_ + ^ f xo T(i)( v )^y 

7T m 2 J " C 2 X + C 2 J Q ±Z {V) y • 



(28) 

(29) 

(30) 
'o y 

Consider a case when the resulting electrical current produced by ex- 
ternal EEF is purely sinusoidal, I z ' = 0. Than it follows from (27), (30) 



c\ 



(r) 



0, £ 



CO 



-<#>. 



(31) 
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At xq <C 1 we neglect logarithmic terms in (27)- (28), and terms with 
integrals in (27)-(30), so the details of the current distribution over the 
cylinder radius are of a little importance. In this approximation 



Cf = ^> £T = ^ln^. (32) 



,( r ) _ TTU {i) {r) _ 2u X (i) 

The solution for the electromagnetic field of the long wave emitted 
by the cylinder with the sinusoidal electric current, starting from the 
surface of the cylinder until the wave zone, follows from (14): 



E z = —z-lk [Yq(x) cos u>t — Jq(x) smut], 
c z 

B^ = ^-Iq [Yi(x)smujt + J\(x) cos art]. (33) 

Near the cylinder we have, with account of the expansions (21) 



E z = — g-ig (— In — coswt — smuit) (34) 

2uj ( j) x 2uj (i) ruj 

~ — =- in In — cos wt = —^-ln In — cos wi, 
c 2 ° 2 c 2 ° 2c 

(i) 

ttlu a) 2 x 2ui a) 2JA 

i><i = ^-jA ( sinwr H — coswi) ~ -=— In smwt = sinwt. 

c J 7TX 2 c 2 x cr 

(35) 
Note, that in the near zone of the long wave the magnetic field adia- 
batically follows the current through the cylinder (Landau and Lifshits, 
1982). The distribution of the electrical field is similar to the one of the 
linearly growing current in the cylinder (Bisnovatyi-Kogan, 2003). At 
large r we obtain from (33), (12) the expanding cylindrical wave 

B^ = -E z = -^-J—I^ S in(x -j -cot) (36) 

c A V txx 4 



-J /W sin -r-ci --. 

c\cr c 4 



6. Electromagnetic energy flux from jet 

Strong electromagnetic wave generated by oscillations may accelerate 
effectively particles at large distances from the nucleus near the jet, as 
well as at larger radiuses (Bisnovatyi-Kogan and Lovelace, 1995). Let 
us estimate the energy flux in the electromagnetic wave radiated by the 
jet of the length I, and radius r$. If n e is the electron density producing 
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Figure 1. Magnetic and electrical fields around the infinite cylinder with the radius 
ro, and low-frequency u -C c/ro, sinusoidal electrical current along the cylinder axis. 
In the near zone electrical and magnetic fields are varying in antiphase, and far from 
the cylinder r S> c/uo the expanding cylindrical electromagnetic wave is formed, with 
E z — —Bcj,. 

the electrical current, with e as electron charge, than, with account of 
(36) the Pointing flux P = ^-[EB] through the cylinder surface is 



2irr lP 



ttIu 
2^ 



J-c\ 



(37) 



For the amplitude of the electrical current along the cylinder radius 
Iq = 7rrgn e ce, we obtain the energy flux from the jet in the form: 



F 



~2 



e lr u>n e 



2-10 4y erg/s- 



/ 



>'o 



lkpc \lp c / 100 yr \10 10 cm 



n r , 



Here T = — is the period of the electromagnetic wave. Part of the 
radiated energy is used for particle acceleration up to very large energies 
(Bisnovatyi-Kogan and Lovelace, 1995), and support the jet radiation 
at different energy regions of electromagnetic spectrum. 



7. Generation of strong electromagnetic wave by proper 
oscillations in a separate blob. 

The mechanism of shock acceleration of particles, often considered 
(Eilek and Hughes, 1990), is not certain, it is unlikely that shock 
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acceleration can give a fairly uniform brightness jet as observed in 
some cases. The mechanisms of magnetic field reconnection (Roman ova 
and Lovelace, 1992) and plasma turbulence acceleration (Eilek and 
Hughes, 1990) are also highly uncertain. When plasma density in the 
surrounding medium is small, the electromagnetic wave generated by 
the nonpotential plasma oscillations of the confined body is emitted 
outside and can accelerate particles. When the emitted wave is strong 
enough it washed out the medium around and the density may become 
very small, consisting only of the accelerated particles. The solution 
of the whole problem of the dynamical behavior of the confined knots 
embedded into the large scale elongated magnetic field and producing 
the toroidal field can be solved by self-consistent calculations of the knot 
oscillations, using together the hydrodynamical and complete Maxwell 
equations. In order to estimate properties of a long wave radiation of 
by oscillating knot we solve instead the idealized problem having the 
analytical solution. Consider linear plasma oscillations of the infinitely 
long uniform cylinder. The problems of such kind have been inten- 
sively studied for plasma wave-guides (Kondratenko, 1976). The main 
difference in this problem is another boundary conditions which suggest 
vacuum state around the cylinder. When considering linear electromag- 
netic oscillations in the static plasma cylinder, only Maxwell equations 
(3) with time dependence in the form ~ exp(— iuit) are needed. The 
background constant field B z = Bq is adopted. Dependence if j on 
E, B and Bq can be obtained, using the expression for the dielectric 
permeability 



47T2 

€ij = Sij H Wij(e) + <7ij(p)], ji = Wij(e) + (Tij(p)]Ej (39) 

Here we consider for simplicity pure hydrogen plasma. The components 
for <Jij{a) in the case of perfect conductivity are (Kondratenko, 1976) 

(Tii (a) = CT22 (a) = ~. n m 2 > ( 40 ) 

1 Wpa^Ba . . % ^ 

(712(a) = -<T 2 i(a) = --; 5 t~ , 0"33(a) = -. 

47T LU Z — 0J Ba 47T ijj 

where uj pa and oJBa ar e plasma and Larmor frequencies of electrons 
(e) and protons (p): uo pa = ^ e ° , uj B a = f^r- I n tne cylindrical 
coordinates (r, 0, z) with ^ = 0, ■£- = ik we have from (3), (39) 

47T XLJ 

-ikB^ = —j r E r , (41) 

c c 
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ikB r 



dB r 
dr 



4tt 



—ikEfk = — B r 



c 
ikE, 



J0 



■lit! 

c 

dE z 
dr 



E r : 



1 d 

r dr 



(rB^) 



4tt 



c 



IUJ 



E z 



iuj „ 1 d , „ . w „ , , 
-% -^(^) = -^- (42) 



We are interested in long wave oscillations with u) <C oj- 
ujb p , so approximately 



pei ^ppi 



^Be 



_ _ 1 Up ] ^ 

Jr ~ 4vr u 2 Bp r ' H 



4vr w| 



Ea>, jz 



J_L0pe 

4-7T a; 



E x . 



(43) 



Substituting (43) into (41), (42) we can see that two types of long waves, 
corresponding to different polarizations exist independently: electric- 
type (E) waves with B z = and nonzero (E r , E z , B^), and magnetic- 
type (B) waves with E z = and nonzero (B r , B z , E^). The equations 
for E— wave has a form 



d 

dr 


- — (rZ E ) 
r dr 




1 d fdE z \ /u 2 e 
r dr\ dr J ' \lo 2 


Ze = E r or B^, and 




c\ = c 2 (l 


+ i") 



c 2 
c 2 



£l*. 



*'* 



0. 



0, 



-i 



^Bp, 



1 + 



4-7T/9C 2 



-1 



(44) 
(45) 

(46) 



is a speed of the Alfven waves. The equations for a B— wave has a form 



d 


\ ld ( 7 


dr 


r dr 


\ d ( dB z 


r 


— r 

dr V dr 



+ 



u 



LO 



+ Ur - k z )B, 



k 2 )Z B = 0, 



0. 



(47) 



(48) 



where Zb = B r or E^, Outside the cylinder equations (44), (45), 
(47), (48) retain their form after substituting c everywhere instead of 
ca- Solutions of these equations are given by Bessel functions: 



Z b ,e = Z BOt EoJi(nB,Er), B z = B Z oJo(K B r), E z = E Z0 J Q (K E r) 
where k% = 4- - k 2 , k\ = { ^- - I 



(49) 
< i i — ? — ' i i —j- — -t i inside the cylinder. 
Outside the same solutions (49) are valid with k 2 = ^ — k 2 , instead of 



—7T- 
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n 2 E or k%. Discreet values of Kb and ke are determined by dispersion 
equation, obtained from boundary conditions on the surface of the 
cylinder. For B— wave the components B r and E^ are continuous on 
the boundary and B z = is taken at the boundary r = tq — 0. That 
leads to the relations 

( A 2 \ 

^B„r = A „ , uJ 2 n = c\ [k 2 + -% . (50) 

For E— wave E z is continuous and E r = at r = ro — 0, so we have 

HEiro = Aii, u;f = c^fc 2 1 + 2 U 2 , (51) 

where Ao re and Ah are zero's of the Bessel functions J (x) and J\{x). 
We have then for the wave vectors outside the cylinder 



-I = 44 - k 2 (i - #) (5), 4 = -k 2 (i - -3d* ) ( ^). 



C 2 Tn V C 2 J ' V i I A ii 



c^ 



'o v ^ / v 1 + ^-,- 

It is clear from (46) that k 2 < for E— waves, but for each k there is n* 
such that k 2 > for n > n* and Alfven 5— wave is transforming into 
electromagnetic one. The vacuum values of B r (ro + 0) and E^r® + 0) 
at the boundary are obtained from the continuity conditions, and the 
outer value of B z (ro + 0) is found from the last equation in (42). The 
jump AB Z = B z (ro + 0) on the boundary determines the induced 
surface current i^. Very long electromagnetic wave emitted with the 
amplitude E^ could accelerate particles up to energies e k> eEr^ pz 
alO 3 ergs for ro = 1 pc and E = a ■ 10~ 6 in CGS. Real energy of the 
accelerated particles could be much less due to radiation losses. 



8. Conclusion 

Similarity of the acceleration mechanisms of the particles in the pulsars 
and relativistic jets could be a reason of the similarity of the high energy 
radiation around 100 Mev which have been observed by EGRET in 
number of radiopulsars (Thompson et al., 1992); quasars and AGN 
(Thompson et al., 1993), but in no other compact objects. 
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